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ABSTRACT 

In this paper we use high quahty X-ray observations from XMM-Newton and Chan- 
dra to gain new insights into the explosion that originated Tycho's supernova 433 years 
ago. We perform a detailed comparison between the ejecta emission from the spatially 
integrated X-ray spectrum of the supernova remnant and current models for Type la 
supernova explosions. We use a grid of synthetic X-ray spectra based on hydrodynamic 
models of the evolution of the supernova remnant and nonequilibrium ionization calcu- 
lations for the state of the shocked plasma. We find that the fundamental properties of 
the X-ray emission in Tycho are well reproduced by a one-dimensional delayed detona- 
tion model with a kinetic energy of ~ 1.2 • 10^^ erg. All the other paradigms for Type 
la explosions that we have tested fail to provide a good approximation to the observed 
ejecta emission, including one-dimensional defiagrations, pulsating delayed detonations 
and sub-Chandrasekhar explosions, as well as deflagration models calculated in three 
dimensions. Our results require that the supernova ejecta retain some degree of chem- 
ical stratiflcation, with Fe-peak elements interior to intermediate mass elements. This 
strongly suggests that a supersonic burning front (i.e., a detonation) must be involved 
at some stage in the physics of Type la supernova explosions. 

Subject headings: hydrodynamics — ISM:individual(SN1572) — nuclear reactions, nu- 
cleosynthesis, abundances, — supernova remnants — supernovae: general — X-rays:ISM 
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1. INTRODUCTION 

The Stella nova of 1572 is the only historical supernova (SN) that can be classified with 
some degree of confidence as Type la based on its observed light curve and color evolution (Ruiz- 
Lapuente 2004). Our knowledge of the universe has improved greatly since this momentous event 
was recorded by Tycho Brahe and other astronomers in the sixteenth century, but now, as then, 
there are many unanswered questions regarding the nature of supernovae. Type la SNe play a key 
role as cosmological probes (Riess et al. 1998; Perlmuttcr ct al. 1999), and they constitute the most 
direct evidence for the accelerating universe (Leibundgut 2001), but despite the continuing efforts 
of theorists during the last decades, our understanding of these objects is far from being complete. 
There is a general agreement that Type la SNe are the result of the thermonuclear explosion of a 
C+0 white dwarf (WD) that is destabilized by accretion from a companion star in a close binary 
system. This so-called 'single degenerate' scenario has been recently substantiated by the discovery 
of a G-type star near the center of the Tycho supernova remnant (SNR) that appears to be the 
surviving companion to the WD that exploded as a SN in 1572 (Ruiz-Lapuente et al. 2004), and by 
the presence of a weak Ha signature in the late time spectrum of the unusual Type la SN 2002ic 
(Hamuy et al. 2003). However, the fundamental details of the explosion itself are still obscure, and 
a number of contending models or paradigms are currently being considered. 

Most of these paradigms assume that the WD becomes unstable and then explodes when its 
mass approaches the Chandrasekhar limit (1.4 Mq). In this case, the properties of the explosion are 
determined mainly by the propagation mode of the burning front through the star, which can be 
supersonic (detonations), subsonic (deflagrations), or a combination of both (delayed detonations, 
pulsating delayed detonations). In the prompt detonation models, the burning front propagates 
supersonically through the star, incinerating almost its entire mass to Fe-peak nuclei. In pure 
deflagration models, the burning front propagates subsonically and the star has the time to expand 
ahead of the flame. When the expansion velocity of the unburnt material becomes comparable to 
that of the burning front, the flame quenches, leaving a mixture of unburnt fuel (C and O) and ashes 
(Fe-peak nuclei) in the ejecta, usually with trace amounts of intermediate mass elements (IMEs: 
Si, S, Ar, Ca, etc.). In one-dimensional deflagration models, the ejecta are layered, with the ashes 
interior to the fuel, but three-dimensional calculations have shown that the turbulent properties of 
the burning front lead to an efficient mixture of fuel and ashes through the ejecta. In the delayed 
detonation models, the burning front begins as a deflagration and then it is artificially forced to 
make the transition to a detonation, usually when a prescribed transition density ptr is reached. 
In this kind of explosions, most of the WD is burnt and very little C and O is left behind, but the 
expansion of the star during the initial deflagration allows for the detonation to burn a sizable mass 
of ejecta at intermediate densities, leading to a significant production of IMEs. Pulsating delayed 
detonations are a variation of this scenario, where the flame quenches after a brief deflagration 
phase and the WD re-collapses, triggering a detonation in the process. The ejecta are stratified 
in all the flavors of delayed detonations simulated in one dimension, with Fe-peak nuclei interior 
to IMEs, and some unburnt C and O in the outermost layers. Whether this is also the case for 
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three-dimensional delayed detonations is still under debate. Given the appropriate circumstances, 

it is also possible for the WD to explode with a sub-Chandrasckhar mass, leading to an entirely 
different class of models whose initial conditions are not so well determined. Sub-Chandrasekhar 
explosions usually involve the ignition under degenerate conditions of a layer of accreted He on 
top of the C+0 WD, which sends a strong pressure wave towards the center of the star that 
eventually triggers a detonation. The structure of the ejecta is complex, with unburnt C and O and 
IMEs sandwiched between two regions rich in Fe-peak nuclei. For more details on these explosion 
paradigms and a complete set of references, see the reviews by Hillebrandt & Niemeyer (2000) 
(mostly onc-dinicnsional models). Bravo et al. (2005) (three-dimensional deflagrations), and the 
recent works by Garcia-Senz & Bravo (2003) and Gamezo et al. (2005) (three-dimensional delayed 
detonations) . 

There have been numerous attempts to constrain these theoretical models, mostly through the 
optical/IR spectra (e.g., Hofiich & Khokhlov 1996; Fisher et al. 1997; Hofiich et al. 1998; Wheeler 
et al. 1998; Branch et al. 2005; Kozma et al. 2005) and light curves (e.g., Hofiich et al. 1998; 
Stritzinger et al. 2005) of Type la SNe, and also through the few available gamma-ray observations 
(Milne et al. 2004, and references therein). Although some preference has been shown for delayed 
detonation models in these studies, it cannot be said that a general consensus has been reached 
in this matter, and the physical mechanism responsible for Type la SN explosions still remains 
an open issue. In the present paper, we address the problem of constraining the physics of Type 
la SN explosions from a different perspective. Instead of focusing on the emission from the SNe 
themselves, we take advantage of the excellent quality of the existing X-ray observations of SNRs 
provided by XMM-Newton and Chandra to probe the physical mechanism responsible for the 
explosions. The necessary theoretical groundwork was laid down in Badenes et al. (2003) and 
Badenes et al. (2005a) (henceforth, Paper I and Paper II), where we showed that the structure 
of the SN ejecta has a profound impact on the thermal X-ray emission from Type la SNRs. As 
a consequence, the X-ray spectra have the potential to pose strong constraints on the kind of 
explosion that originated the SNR. Here we apply the models developed in Papers I and II to the 
X-ray spectrum of the Tycho SNR. We begin by summarizing the observations of this object in § 2. 
Based on the observational properties of Tycho, we outline our strategy for comparing our models 
to the observations in § 3. In § 4, we extract spatially integrated spectra from large regions of 
an XMM-Newton observation of Tycho and derive the fundamental properties of the line emission 
from the shocked SN ejecta. We compare our models with these observations in two stages: first 
we select the most promising models based on the properties of their line emission in § 5 and then 
we compare the best candidates to the spatially integrated spectrum in § 6. In § 7 we consider 
the spatial distribution of the line emission from a Chandra ACIS-I observation, and examine the 
implications for our models. Finally, we discuss the performance of the models from a global point 
of view in § 8, and we summarize our conclusions in § 9. 
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2. OBSERVATIONS 

2.1. The supernova of 1572 

As early as 1945, Baadc used the observations made and compiled by Tycho Brahe to derive 
a light curve that enabled him to identify the stella nova of 1572 as a Type I supernova (Baade 
1945). Since then, the classification of the supernova as determined from its light curve has been 
controversial, and some authors have claimed that SN 1572 was a subluminous event (see van den 
Bergh 1993; Schaefer 1996, and references therein). Ruiz-Lapuente (2004) performed a detailed 
reanalysis of the sixteenth century records and evaluated the uncertainties in the data in order to 
describe the light curve in the terms used nowadays to characterize Type la SNe. According to 
Ruiz-Lapuente, SN 1572 was a normal Type la SN within the uncertainties associated with the 
data (which are large), with a peak visual magnitude of My = —19.24 — 5log{D / 3. Okpc) it 0.42 
mag, and a stretch factor ^ s ^ 0.9. The author found an extinction of Ay = 1.86 it 0.12 mag 
and an average reddening of E{B — V) = 0.6 it 0.04 in that direction of the sky, which combined 
with the value of Mv yield an estimate of 2.8 it 0.4 kpc for the distance to the SN, assuming 
i7o = 65km-s-i-Mpc-^ 



2.2. The Tycho SNR: Radio and optical observations 

In the radio continuum observations at 1375 MHz (A = 21 cm), the Tycho SNR appears as a 
clearly defined shell with an approximate angular diameter of 8'. The shell is nearly spherical from 
the northwest to the southeast, with an irregular outbreak accompanied by a slight brightening 
to the north, northeast and east. Reynoso et al. (1997), used VLA observations to estimate the 
expansion parameter of the forward shock, rifyj^ (defined as rj = dlog{r fyj^) / dlog{t)) , along the rim 
of the SNR shell. They found an average of r)fyjd = 0.47 ±0.03, with distinctly lower values towards 
the north and east. It was found later that the lower value of r//.u,'(i in the northeast is due to an 
interaction with a dense cloud of neutral hydrogen (Reynoso et al. 1999; Lee et al. 2004). There 
is a puzzling, and still unsolved, disagreement between these radio measurements and the value of 
Vfwd = 0-71 ± 0.06 inferred from X-ray observations by Hughes (2000). The expansion parameters 
of several interior features, measured using the same techniques, are consistent in radio {rf c± 0.44, 
Reynoso et al. 1997) and X-rays {rj ~ 0.45, Hughes 2000). 

At optical wavelengths, only a few faint filaments of Balmer line emission from H are visible 
at the outer rim of the SNR. No evidence for any optical emission other than the Balmer H lines 
was found in the interior of the SNR (Kirshner & Chevalier 1978; Smith et al. 1991; Ghavamian 



^The stretch factor s is a parameterization of the hght curve width/shape relationship, defined as the linear 
broadening or narrowing of the rest-frame time scale of an average template light curve that is required to match 
the observed light curve. Values of s close to 1 denote 'normal' Type la SNe. See Perlmutter et al. (1997) for more 
details. 
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et al. 2000). In particular, optical emission from radiatively cooled ejecta has not been detected. In 
Ghavamian et al. (2001), the spectrum of the brightest knot in the eastern rim (knot g from Kamper 
& van den Bcrgh 1978) was examined in detail, yielding a velocity between 1940 and 2300km • s~^ 
for the forward shock. This bright knot, however, is situated in the eastern rim, where the SNR 
is interacting with denser material, and its properties might not be representative of the overall 
dynamics of the blast wave. 

As is often the case with Galactic SNRs, the distance to Tycho is very uncertain. Different 
techniques yield different and even contradictory results ranging anywhere between 1.5 and 4.5 kpc, 
with most estimates converging around 2.5 kpc (see Schaefer 1996, for a review). Among these, we 
believe that the most reliable results are obtained through the combination of velocity estimates 
from knots in the Balmer-dominated forward shock with optical proper motions measured over 
long temporal baselines. We will adopt the range of 1.5-3.1 kpc obtained by Smith et al. (1991) 
with this method as a conservative estimate for the distance, and we note that the value of 2.8 kpc 
found by Ruiz-Lapuente (2004) for SN 1572 falls within this range. 

2.3. The Tycho SNR: X-ray observations 

The most fundamental properties of the X-ray emission from the Tycho SNR have been known 
since the time of the ASCA satellite (Hwang & Gotthelf 1997; Hwang et al. 1998), but the quality 
of the available data has increased dramatically thanks to XMM-Newton (Decourchelle et al. 2001) 
and Chandra (Hwang et al. 2002; Warren et al. 2005). Morphologically, the X-ray emission from 
the Tycho SNR has a shell-like structure with a thin rim of high energy continuum emission. This 
rim is coincident in position, but not in brightness, with the outer edge of the radio emission, and 
traces the position of the forward shock. The regions interior to this rim show very strong emission 
lines of Si, S, Ca, Ar and Fe, which are thought to arise from the shocked supernova ejecta. 

Several significant results about the dynamics and X-ray emission of the Tycho SNR were 
presented in Warren et al. (2005), where a principal component analysis technique was applied to a 

~ 150 ks Chandra ACIS-I observation. Warren et al. determined the positions of the forward shock 
(FS), reverse shock (RS) and contact discontinuity (CD) between ejecta and AM, and derived ratios 
for their average radii of 1:0.93:0.71 (FS:CD:RS). Because the CD is so close to the FS, these ratios 
are incompatible with one-dimensional adiabatic hydrodynamics with 7 = 5/3 for the shocked AM. 
Three explanations for the presence of ejecta close to the FS can be found in the literature: the 
presence of clumps with a high density contrast in the ejecta (Wang & Chevalier 2001), the effect 
of cosmic ray (CR) acceleration at the FS (Decourchelle et al. 2000; Blondin &: Ellison 2001), and 
the interaction of Rayleigh- Taylor fingers with circumstellar cloudlets (Jun et al. 1996). The last 
mechanism is unlikely to be relevant in the case of the Tycho SNR, because there is no known process 
that would produce an almost isotropic distribution of cloudlets in the circumstellar environment 
of a Type la SN progenitor (see Figure 4 in Warren et al. 2005, for the azimuthal distribution of 
FS and CD radii). According to Warren et al., the absence of prominent variations in the spectral 
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properties of the ejecta along the CD seems to favor CR acceleration over ejecta clumping as the 
dominating mechanism. 

Warren et al. also found that the properties of the X-ray emission from the shocked AM are 
compatible with a CR-dominated blast wave. On one hand, the emitting region behind the FS 
is much thinner than would be predicted by the thermal emission behind an adiabatic shock (see 
also Ballet 2005). On the other hand, the featureless spectra extracted from the shocked ambient 
medium (AM) were well fitted by a power law with an index of ~ 2.7. This value is consistent 
with the index of 2.72 found in the high energy Ginga observations by Fink et al. (1994). Fits 
with thermal models at fcT ~ 2 keV are also possible, but they require extremely low values for the 
ionization timescale (a few times 10^ cm~^ • s, Hwang et al. 2002), which are incompatible with the 
basic properties of the Tycho SNR (see the discussion in § 7.4.1 of Warren et al. 2005) 

The emission from the shocked ejecta, as seen in the spectral bands corresponding to the most 
prominent emission lines from heavy elements, has a shell-like morphology, with the Fe Ka line 
image appearing more diffuse and peaking at a smaller radius than the others (Hwang &; Gotthelf 
1997; Decourchelle et al. 2001). The Fe Ka emission is not only spatially distinct from the other lines 
(including the Fe L complex); it also has different spectral properties, with a higher temperature 
and a lower ionization timescale (Hwang et al. 1998). The apparent symmetry of the X-ray line 
emission and the absence of significant Doppler shifts suggests an overall spherical geometry, but 
local inhomogeneities in the line emission are manifest, particularly in the form of bright clumps 
in the southeast (see Vancura et al. 1995; Decourchelle et al. 2001). 

3. GOALS AND STRATEGY 

The goal of this paper is to model the thermal X-ray emission from the shocked ejecta in the 
Tycho SNR using the grid of synthetic spectra presented in Papers I and H. Our objective is to use 
these synthetic spectra to constrain the physics of the Type la SN explosion that gave birth to the 
remnant in 1572. We shall make no attempt to model the featureless emission from the shocked 
AM. Instead, we adopt the view expressed in Warren et al. (2005) that this emission is predom- 
inantly nonthermal, and that the dynamics of the FS are strongly modified by CR acceleration. 
Although the one-dimensional hydrodynamic calculations with 7 = 5/3 that underlie our synthetic 
spectra cannot reproduce the properties of the shocked AM or the dynamics of the FS under these 
conditions, this is not a serious drawback. Several fundamental quantities like the temperature of 
the shocked AM or the expansion parameter of the FS arc not known with a significant degree 
of accuracy (see the discussions in § 2.2 and § 2.3). Rather than use these poorly determined 
quantities to constrain the hydrodynamic models underlying our synthetic spectra, we attempt to 
reproduce the observed X-ray emission first and then verify that the SNR dynamics is compatible 
with the known properties of the AM and FS. 

Regardless of the conditions at the FS, we believe that our hydrodynamic models with 7 = 5/3 
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are a good approximation for the dynamics of the shocked ejecta. Diffusive shock acceleration is 
mediated by strong magnetic fields, and leads to high compression ratios and low temperatures 
in the shocked plasma (see Ellison et al. 2004; Ellison et al. 2005). From an observational point 
of view, both the position of the RS measured by Warren et al. (2005), deep into the ejecta and 
away from the CD, and the strong Fe Ka emission that originates from the hot plasma just behind 
the RS are strong arguments against efficient CR acceleration at the RS. Prom a theoretical point 
of view, it is hard to see how the necessary magnetic fields could have survived for hundreds of 
years inside the freely expanding SN ejecta, even making the most optimistic assumptions about 
the initial magnetic field strength in the progenitor systems of Type la SNe. 

As we discussed in Paper II (§ A. 2), a detailed comparison between our synthetic spectra and 
observations is not straightforward. On one hand, the uncertainties in the atomic data and other 
factors do not make the model spectra amenable to the usual fitting techniques, and on the 
other hand, the parameter space is very large even for one-dimensional calculations. In the present 
work, we solve these problems by taking advantage of the characteristics of the X-ray emission from 
the Tycho SNR to devise an efficient comparison strategy. First, we extract spatially integrated 
spectra from large regions in the SNR, and determine the fundamental properties of the high- 
energy (> 1.5 keV) line emission (§ 4). As we have seen in § 2.3, this line emission comes entirely 
from the shocked ejecta, so we use the observed line flux ratios and line centroids to reduce the 
dimensionality of the problem. We require that our synthetic X-ray spectra reproduce as many of 
these flux ratios and centroids as possible, and in the process we identify the most successful Type 
la SN explosion models and the most promising regions of the associated parameter space (§ 5). 
Finally, we compare the best models to the entire spatially integrated spectrum, applying several 
consistency checks to the requirements of each model, such as interstellar absorption, AM flux, and 
derived distance estimates (§ 6). 



4. SPATIALLY INTEGRATED X-RAY SPECTRUM AND LINE EMISSION 

Among the publicly available X-ray observations of the Tycho SNR, the best spatially inte- 
grated spectra are provided by the EPIC MOS cameras on board XMM-Newton. The observation 
analyzed in Decourchelle et al. (2001) has a total exposure time of about 12 ks for each of the 
EPIC MOS CCDs, providing spectra of good quality for this bright object ^. Since our models are 
one-dimensional, we focus on the western sector of Tycho, avoiding the prominent inhomogeneities 
in the ejecta emission found in the southeast and the interaction with a denser AM towards the 
northeast (see discussion in § 2.2 and § 2.3). In the west, the ejecta emission is homogeneous, and 
the FS has a smooth, almost circular shape. Using the standard XMM-Newton science analysis 
system, we have selected two spectral extraction regions in the EPIC MOSl image (see Fig. 1). 



^Tho data arc available from the XMM-Newton science archive (XSA): 
http: / /xmm. vilspa.esa.es / external/xmm_data_acc/xsa/ . 
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Region A covers the entire western sector between two major outbreaks in the FS in the north and 

south; it corresponds loosely to region V in Reynoso et al. (1997), with a range in position angle 
of 200 < < 345 degrees (counterclockwise from the north). Region B is a subset of region A, 
covering the range 200 < 6 < 300, and has been selected to exclude an area in the northwest that 
shows an increased flux in the Fe and Si emission (see Figs. 2 and 3 in Decourchelle et al. 2001). 
The spectra extracted from regions A and B are shown in Figure 2. 

Many of the lines and line complexes in the spectra of Figure 2 are blended due to the limited 
spectral resolution of CCD cameras in the X-ray band. In order to characterize the properties of 
the line emission in regions A and B, we have fitted the extracted spectra between 1.65 and 9.0 
keV with a model adapted from Hwang &; Gotthelf (1997), consisting of fourteen Gaiissian lines 
plus two components for the continuum (a thermal bremsstrahlung and a power law), attenuated 
by the interstellar absorption. For the neutral hydrogen column density, we have adopted a value 
of Nff = 0.6 • 10^^ cm~^, which is compatible with the 0.55 — 0.59 • 10^^ cm~^ range found by Hwang 
et al. (1998) from fits to the integrated ASCA spectrum. It is also close to the values determined 
by fitting the Chandra spectra in several locations along the western rim with nonthermal models: 
0.58-0.69-1022 cm-2 (Warren et al. 2005); 0.71-0.79-1022 cm~2 (Hwang et al. 2002). Conversion 
of the value of Ay obtained by Ruiz-Lapuentc (2004) to Nh with the Ay-Nfj relation of Predehl 
&; Schmitt (1995) yields a lower column density of 0.34- 1022 cm~2^ but the uncertainties associated 
with this estimate are large. We stress that the interstellar absorption in the X-ray band is hard 
to constrain in objects with complex spectra such as SNRs, and variation of Nh across the surface 
of an object as extended as Tycho cannot be discarded (sec § 5.2 for a discussion). For the power 
law index, we have adopted the value of 2.72 found by Fink et al. (1994) and confirmed by Warren 
et al. (2005). 

The fits to regions A and B are shown in Figure 3. The statistics for the fit to region A are 
^2 = 462.04 with xV^OF = 1.52, and for region B = 373.75 with xV^OF = 1.38. The 
fitted temperatures for the thermal bremsstrahlung component are 0.48 and 0.40 keV, respectively. 
We note that these values are rather different from those obtained by Hwang & Gotthelf (1997) 
assuming a purely thermal continuum (these authors used a main component at 0.99 keV, with 
another at 10.0 keV to reproduce the high energy continuum). The fourteen lines included in the 
model, along with the fitted centroids and fluxes, are given in Table 1, where the common notation 
of a, P, and 7 has been used to label the lines corresponding to transitions from quantum levels 
2, 3, and 4 to level 1. The quality of the data set did not allow the centroids of the weakest lines 
to be fitted independently, so these parameters were fixed. The He7/He/3 line flux ratios of Si, S 
and Ar, and the Si Ly/3/Si Lya ratio, were also flxed in the fits to the values listed in Table 1, 
which correspond to the values at T = lO''' K This allows for an adequate (if simple) treatment 
of these blended lines, and is justified because the fiux ratios vary by only 10%-20% over a decade 



^Some of these flux ratios have been updated with respect to the values listed in Hwang & Gotthelf (1997) using 
the atomic data from the ATOMDB compilation available at http://cxc.harvard.edu/atomdb/. 
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in temperature (for details, see Hwang &; Gotthelf 1997). The most important line flux ratios are 
listed in Table 2. For these flux ratios, we have used the S Hea blend as a reference instead of Si 
Hea in order to make their values less dependent on A^^^. 

An inspection of Tables 1 and 2 reveals that the properties of the line emission in the western 
half of Tycho are remarkably uniform. Indeed, the only noticeable differences between the two 
extraction regions are in the centroids of the Ar Hea, Ca Hea and Fe Ka blends, which are a few 
eV lower in region B (the deviation is below 0.5% in all cases, and within the 90% confldence ranges 
for Ca Hea and Fe Ka). This suggests a lower ionization state for these elements in region B, but 
the origin of this lower ionization and its relationship to the brightening in the ejecta emission 
in the northern part of region A is unclear. In the absence of a plausible physical interpretation 
for these effects, we use region B as representative of the shocked ejecta emission in the reminder 
of the paper, but we stress that the differences between regions A and B are very small. Our 
XMM-Newton results for the western half of Tycho are very similar to those obtained by Hwang 
&; Gotthelf (1997) for the ASCA observation of the entire SNR. Although the measurements were 
made in different regions of the SNR, assuming different values for Nh, with different models for 
the continuum, and with different instruments, the error bars overlap in almost all the line flux 
ratios and centroids, and the deviations arc never larger than a few percent. The only case where 
this is not true is the centroid of the Ca Hea blend (3.879 and 3.867 keV in regions A and B vs. 
3.818 keV in Hwang & Gotthelf 1997). The neighboring Ar He/3 line is stronger in our fits, shifting 
the Ca Hea centroid towards higher energies that are closer to the expected value of 3.88 keV. 



5. MODELING THE LINE EMISSION 

5.1. Synthetic spectra 

The method to generate the synthetic X-ray spectra that we use to model the line emission 
from the Tycho SNR is described in Papers I and II. For each Type la SN explosion model, we 
simulate the interaction of the ejecta with a uniform AM using a one-dimensional hydrodynamic 
code. We assume that the plasma is a nonrelativistic monoatomic ideal gas with 7 = 5/3. Then 
we calculate the nonequilibrium ionization (NEI) and associated processes for each fluid element in 
the shocked ejecta, taking as inputs the temporal evolution of density and specific internal energy 
from the hydrodynamic simulation and the chemical composition from the SN explosion model. 
To perform these calculations, we use a method adapted from Hamilton & Sarazin (1984), which 
takes into account important details that are relevant to the specific case of a plasma dominated 
by heavy elements. In particular, the influence of the large number of electrons ejected during the 
ongoing ionization on variables like the plasma temperature (noted by Brinkmann et al. 1989) is 
dealt with in a self-consistent way. Once the electron temperatures and ionization states in all the 
shocked ejecta are calculated in this way, a spatially integrated X-ray spectrum is produced using 
an appropriate spectral code. 
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This simulation scheme does not take into account radiative or ionization losses or thermal 
conduction in the shocked plasma. Hamilton & Sarazin (1984) showed that ionization and radiative 
losses arc of the same order in a plasma dominated by heavy elements. In order to assess the 
importance of these processes, we applied the post-facto estimation of radiative losses described 
in section 3.5 of Paper I and section A. 3 of Paper II to all the models presented here. We found 
that radiative losses are unimportant for an object of the age of the Tycho SNR, with only a 
few inconsequential exceptions that we will specify in § 5.3. This is in agreement with both the 
conspicuous absence of optically emitting ejecta (see § 2.2) and the results of simulations that do 
include radiative and ionization losses, like Sorokina et al. (2004). These authors also showed that 
efficient thermal conduction destroys any temperature gradients on timescales that are shorter than 
the age of the Tycho SNR. Thus, the presence of efficient thermal conduction is incompatible with 
the spatial morphology of the Fe Ka and Fe L emission in the Tycho SNR (see § 2.3), which implies 
that a temperature gradient must exist in the shocked ejecta. 

The synthetic spectra for the shocked SN ejecta are characterized by four parameters only: 
the Type la SN explosion model used, the age of the SNR t, the density of the AM pam, and the 
ratio between the specific internal energies of the electrons and ions at the reverse shock, /? = Sej^i, 
which represents the efficiency of coUisionless electron heating. The dependence of the ejecta X-ray 
emission on each of these parameters is discussed in Papers I and II. 

The 32 Type la SN explosion models presented in Papers I and II include examples of all the 
paradigms described in § 1. Among these, we have selected a more convenient subgrid of 12 models: 
three delayed detonations (DDTa, DDTc, and DDTe), three pulsating delayed detonations (PDDa, 
PDDc, and PDDe), three deflagrations (DEFa, DEFc, and DEFf), one prompt detonation (DET), 
one sub-Chandrasekhar explosion (SCH), and one deflagration calculated in three dimensions with 
a smooth particle hydrodynamics code (B30U from Garcfa-Senz & Bravo 2005). In order to bench- 
mark the results obtained using this subgrid of models, we have chosen to add three well-known 
Type la SN explosion models calculated by other groups: model W7 (Nomoto et al. 1984), a 'fast 
deflagration' that has become a widespread standard for comparison with observations of Type la 
SNe, model 5p0z22.25 (Hofiich et al. 2002), a delayed detonation calculated with a resolution 4 
times greater than the DDT models of our grid, and model b30_3d_768 (Travaglio et al. 2004), a 
deflagration calculated in three dimensions with an Eulerian code and a high spatial resolution. 

Of the three remaining parameters in the synthetic spectra, the age t is known to be 433 years 
for the Tycho SNR (428 years at the time of the XMM-Newton observation). Thus, only pam and 
(3 can be varied for each of our Type la SN explosion models. We have sampled this parameter 
space with five points in pam ( 2 • lO'^^ , 5 • lO'^^, lO'^^, 2 • lO'^^, and 5 • lO'^^ g • cm'^ ) and 
three points in (3 [fimin^j 0.01, and 0.1), setting t to 430 yr in all the models. These ranges are 
selected to encompass the likely values of (3 (see discussion in § 2.2 of Paper II) and pam in the 



^Defined as Pmin ~ Zame/rrii, with Zs the mean preshock ionization state, me the electron mass and mi the 
average ion mass, see § 2.2 in paper II. 
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Tycho SNR. For each Type la SN explosion model, 15 synthetic SNR spectra have been generated, 
for a total of 225 synthetic spectra. 

5.2. Comparing the models to the observed line emission 

We have selected thirteen observable parameters as diagnostic quantities for the comparison 
between our models and the observations: the eight line flux ratios listed in Table 3 (He(/3 + 7)/Hea 
and Lya/Hea ratios for Si and S, and the ratios of the Hea blends of Si, Ar, and Ca, and the Ka 
blend of Fe, relative to S Hea), plus the centroids of the Si Hea, S Hea, Ar Hea, Ca Hea, and Fe Ka 
blends. In the synthetic spectra, the properties of the line emission are sensitive to variations in 
PAM ; and the Type la SN explosion model. The model with 14 Gaussian lines used in § 4 to derive 
the observed values for the thirteen diagnostic quantities cannot be applied to these synthetic 
spectra, because some of its underlying assumptions (like the fixed value for the He7/He/3 flux 
ratios) might break down in extreme cases. However, the calculation of line fluxes and centroids in 
the synthetic spectra is straightforward if it is performed before convolution with an instrumental 
response. In this format, the lines that contribute to a given blend can be singled out in most cases 
without the risk of contamination from neighboring lines, and the continuum can be subtracted 
easily. The selection energy windows for each of the lines in the unconvolved model spectra are 
listed in Table 3. 

We have verified that these two methods used to determine line fluxes and ratios in the observed 
spectra and in the synthetic spectra give consistent results. To do so, we have generated fake 
data from one particular model with line emission similar to that of the Tycho SNR (from those 
selected in § 5.3), and we have fitted the fake data with the 14 Gaussian line model. Since our 
synthetic spectra only include thermal emission, we substituted the power law continuum by thermal 
bremsstrahhmg in the fit. The line centroids derived from the Gaussian fit were indistinguishable 
from those obtained directly from the unconvolved synthetic spectrum, except for Ar Hea, where 
a ~ 0.4% deviation was observed. The deviation in the line fluxes never exceeded 7%, except again 
in the case of Ar Hea, where it was 11%. 

Before we proceed with the comparison between models and observations, it is important to 
review all the possible sources of uncertainty in the line fluxes and centroids, both for the observed 
spectrum and for the synthetic spectra: 

Instrumental Limitations Several issues related to instrument calibration and data processing 
can limit the accuracy of the energies and fluxes obtained from X-ray spectra, but these limits are 
often hard to estimate. For the instruments on board XMM-Newton, Cassam-Chenai et al. (2004) 
list the deviations found in the centroids of strong lines in a single observation of the Kepler SNR. 
The maximum deviations between the values measured by the difl'erent EPIC cameras (MOSl, 
M0S2, and PN) were of 13.5 eV (0.7%) for Si Hea, 16.1 eV (0.7%) for S Hea and 28.1 eV (0.4%) 
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for Fe Ka (see Table 2 in Cassam-Chenai et al. 2004). Cassam-Chenai et al. also give a value of 
~ 10% for the flux loss in single events, which should be comparable to the value in the Tycho SNR. 
We use these deviations as estimates of the instrumental accuracy for the XMM-Newton spectra. 

Interstellar Absorption The unabsorbcd line fluxes determined from the observed spectra de- 
pend on the adopted value of the interstellar absorption. Varying Nh between 0.4 and 0.8-10^^ cm~^ 
leads to 10% deviations in the Si Hea flux with respect to the value obtained with our fiducial 
hydrogen column density of 0.6 • lO^^cm"^. For the S Hea blend, the deviations are ~ 5%, and 
they become smaller for lines above 3 keV. The Si Hea/S Hea flux ratio, however, is never affected 
by more than ~ 5%, because the flux variations are always correlated. 

Line Extraction and Overlap Due to the complexity of the line emission in the thermal spectra 
of SNRs, it can be difficult to isolate the contribution of a single line. In the observed spectrum, 
the Ar Hea and Ca Hea blends are most affected by this due to the presence of neighboring weak 
lines (e.g., S He7 and Ar He7) that cannot be measured directly^. In the synthetic spectra, this 
problem can be minimized by using the extraction windows listed in Table 3. The only conflict 
that cannot be resolved in this way is in the Ca Hea blend, which completely overlaps with Ar 
He7. In any case, this should not have a major impact on the derived Ca Hea fluxes and centroids 
for the synthetic spectra, because the flux in the Ar Hc7 line only amounts to a few percent of the 
total flux in the Ca Hea energy window (see Table 1). 

Atomic Data In Papers I and II, we used the Hamilton & Sarazin (HS) code to calculate synthetic 
spectra (for a description of the original code and the most important updates, see Hamilton et al. 
1983; Borkowski et al. 2001). This code is implemented in XSPEC, and several simple NEI spectral 
models in this software package make use of its atomic data by default (this is known as version 
1.1 of the NEI atomic data, or NEI vl.l for short). However, the HS code is not adequate to 
model the high quality XMM-Newton spectrum of Tycho, and we used the atomic data from 
version 2.0 of the NEI models in XSPEC instead. These atomic data arc from the ATOMDB data 
base available at http://cxc.harvard.edu/atomdb/ (for a description, see Smith et al. 2001). We 
augmented this atomic data base by adding inner-shell processes, which are missing in NEI v2.0, 
but are important for transient plasmas. We used the most recent published atomic calculations 
for inner-shell processes, such as the K-shell electron impact excitation rates or the radiative and 
Auger rates for K-vacancy states for Fe^^^ to Fc"''^^ from Bautista et al. (2004) and Palmcri ct al. 
(2003). Because the relevant ionization and excitation cross sections for many ions of interest are 
not available, extrapolation along electronic isosequences was often necessary. The quality of inner- 
shell atomic data varies greatly among ions, resulting in potentially large errors that are hard to 



■'This might be the origin of the small discrepancies in the centroid and flux of Ar Hea mentioned above between 
the Gaussian hue fit to the fake data and the windowing of the unconvolved synthetic spectrum. 
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estimate at this time. There are also other problems with the atomic data from ATOMDB in NEI 
v2.0. This data base was originally designed for plasmas in collisional ionization equilibrium, and 
sometimes line cmissivities are not reliable under extreme NEI conditions. For example, we find 
significant differences in the L-shell emission from Fe between NEI vl.l and v2.0 under conditions 
relevant for the shocked ejecta in the Tycho SNR, although both codes are based on the atomic 
calculations of Liedahl et al. (1995) and should produce similar results. This is of no consequence 
for the high energy line emission discussed here, but it will become important when the entire 
spectrum is considered in § 6. In the particular case of Fe L, we use the atomic data from NEI 
vl.l. 

Doppler EfTect The Doppler effect has not been taken into account in the generation of the 
synthetic spectra. Significant Doppler shifts affecting the entire SNR are not expected from the 
low bulk velocities (> —7 ■ 10^ cm • s~^) in the receding environment of Tycho determined by Lee 
et al. (2004) . Regarding shifts associated with the dynamics of the SNR, they should be minimal 
in our spatially integrated spectra, which cover very large regions of the SNR. It is possible that 
some of the lines in the XMM-Newton spectrum are Doppler broadened to some extent (see § 8.3), 
but this kind of effect is very hard to study in CCD data. We shall not pursue this line of research 
in the present work. 

From the preceding discussions, it is clear that a perfect agreement between models and ob- 
servations is not to be expected. In general, we have found that the statistical and instrumental 
uncertainties that affect the observational data are rather small (a few times 10% in the fluxes and 
below 1% in the energies), but the systematic uncertainties that affect the models can be much 
larger, and in many cases they are impossible to quantify. In this context, the tolerance ranges on 
the observed diagnostic quantities should be at once generous enough to provide some flexibility 
in the comparisons and restrictive enough to discriminate the most successful models. For the 
line flux ratios, we have settled on a factor two range above and below the observed value. An 
exception was made for the Si Lya/Si Hea and S Lya/S Hea ratios, where only upper limits can 
be confidently derived from the observations. For the line centroids, we imposed a tolerance range 
of 2% around the observed value. In the case of Fe Ka, which is an isolated line with well known 
atomic data and has the centroid at the highest energy, we imposed a more stringent requirement 
of 1% around the observed value. These tolerance ranges are rather crude, but they are sufficent 
for our present goal of distinguishing between different explosion models. 

5.3. Results for the line emission 

It is not necessary to present here a complete comparison between the values of the diagnostic 
line flux ratios and centroids in the 225 synthetic ejecta spectra and the observed values. Instead, we 
concentrate on six Type la SN explosion models (DDTc, PDDa, DEFc, W7, SCH, and b30_3d_768) 
that are representative of all the pertinent features, and we display their behavior across the pam, 
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(3 parameter space in Figures 4, 5, and 6. In these plots, it is easy to identify which models are able 

to reproduce most of the diagnostic quantities within the specified tolerance ranges (areas shaded 
in gray) for a given combination of pAM and [3. We have repeated this process with each of the 15 
Type la SN explosion models listed in § 5.1, identifying the best values of pam and (3 and noting 
how many of the diagnostic quantities are reproduced - the results are given in Tables 4 and 5. 

Before discussing each class of models in detail, we make a few general comments about the 
results seen in Figures 4, 5, and 6. First, the choice of values for pam and (3 seems adequate 
in that reasonable results are obtained in the middle range, if at all, while the extrema can be 
discarded in most cases. Variations in pam affect the line emission of all elements in most models, 
while variations in f3 affect primarily the Fe Ka emission (see discussion in § 2.4 of Paper II). As 
PAM increases, so does the emission measure averaged ionization timescale {net) in the plasma. 
This moves most line centroids towards higher energies, and increases line fluxes. In particular, 
the Lycc/Hea ratios of Si and S pose very important constraints on the maximum allowed pam- 
Further constraints come from the Ca Hea centroid, which is displaced towards lower energies at 
low PAM in most models, and from the Fe Ka centroid, which is displaced towards Fe Hea (~ 6.7 
keV) at high pam in most models. An increase of [3 leads to higher electron temperatures towards 
the reverse shock (see § 2.2 and § 2.4 in Paper II). This enhances the cmissivity of the Fe Ka 
line, but also inhibits collisional ionization, displacing the Fe Ka centroid towards lower energies. 
Within these general trends, the detailed behavior of each class of models has its own peculiarities: 

Delayed detonations (DDT) This is by far the most successful class of models. In the case 
of model DDTc (see Fig. 4), 12 of the 13 diagnostic quantities are reproduced for pam = 2.0 • 

10~2^ g • cm~^ and (3 between 0.01 and 0.1. The only discrepant quantity is the centroid of the Ca 
Hea blend, whose energy is slightly underpredicted. The other DDT models from the grid, DDTa 
and DDTc, have a very similar behavior, as shown in Tables 4 and 5. The off-grid delayed detonation 
5p0z22.25 is also very successful, but it requires a higher AM density of 5.0 • 10~^^ g • cm~^. 

Pulsating delayed detonations (PDD) This class of models is also relatively successful, spe- 
cially models PDDa (sec Figure 4) and PDDc, which manage to reproduce 12 and 11 of the 13 
diagnostic quantities, respectively. As was discussed in § 4.2 of Paper I, PDD models have chemi- 
cal composition profiles that are very similar to DDT models, but their density profiles are much 
steeper. This leads to higher values of (ngt) in the shocked plasma so that the Lya/Hea ratios of 
Si and S become incompatible with the observations for pAM above 2.0 • 10^^^ g • cm^^. The high 
ionization timescales in the ejecta also affect the Fe Ka centroid, which moves rapidly to higher 
energies as Fe Hea lines become more prominent. We note that PDD models provide a better 
approximation to the Ca Hea centroid than DDT models. 

One-dimensional deflagrations (DEF) The one-dimensional deflagrations from our model 
grid have several serious shortcomings. As can be seen in Figure 5 for model DEFc, many important 
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lines only appear at the highest values of pam-, and even then the properties of the Si and S line 
emission are hard to reconcile with those of other elements. Models DEFa and DEFf behave 
similarly to model DEFc, and none of them can reproduce more than 9 of the 13 diagnostic 
quantities at once, even under the most favorable conditions. Models DEFc and DEFf at pam = 
5.0 • 10~^^ g ■ cm""^ have the peculiarity of being the only models where radiative losses are of some 
importance at the age of Tycho. In these two models, approximately 0.05 M© of C and O in the 
outermost ejecta might undergo runaway cooling. This is only a minor effect, and will not modify 
our conclusions, because the line emission from the other elements is not affected. 

One-dimensional 'fast' deflagrations (W7) It is interesting to analyze the performance of 
this model, because it has been used in attempts to reproduce the X-ray emission from the Tycho 
SNR by several authors (Hamilton et al. 1986; Itoh et al. 1988; Brinkmann et al. 1989; Sorokina 
et al. 2004) . All these previous works used NEI calculations coupled to hydrodynamic models with 
varying degrees of sophistication, and all of them came to the same conclusion: the W7 model does 
not have enough Fe in the outer layers of ejecta to reproduce the Fe Ka flux in the Tycho SNR. 
Our simulations without collisionless electron heating {[3 = (3min) are in agreement with this, as 
shown by the absence of diamonds in the corresponding panels of Fig. 5. Nevertheless, we note 
that partial collisionless heating of electrons at the reverse shock (/3 = 0.1) can explain the Fe Ka/S 
Heo; flux ratio for a wide range of pam values (see the squares in the same panels of Fig. 5). In 
any case, model W7 can only reproduce 9 of the 13 diagnostic quantities, and must be discarded. 

Sub-Chandrasekhar explosions (SCH) Sub-Chandrasekhar explosions lead to a complex 
structure in the ejecta, with an outer region dominated by Fe-peak nuclei bound between two 

strong density peaks (see § 2 in Paper I). This results in high values of (rigt), which in turn restrict 
the range of pam that can reproduce both the Fe Ka centroid and the Lya/Hea ratio of Si and S 
(see Fig. 6). At these low ambient medium densities, the emission measure of Ar and Ca in the 
shocked ejecta is not high enough to account for the observed line emission. 

Prompt detonations (DET) These models are known to be unrealistic for Type la SNe, be- 
cause the ejecta are essentially pure Fe with only trace amounts of intermediate mass elements. 
Not surprisingly, they fail to reproduce the line emission from all the elements in the SNR spectrum 
except Fe. We do not show plots comparing model and observations in this case. 

Three-dimensional deflagrations (B30U, b30_3d_768) These models represent the most 
sophisticated self-consistent calculations of Type la SN explosions that have been performed to 
date. Model b30_3d_768 (see Figure 6) is more successful in reproducing the line emission from 
Tycho than model B30U due to its higher Ar and Ca content. Yet, its predictions are clearly in 
conflict with the observations, mainly due to the presence of large amounts of Fe in the outer ejecta, 
where the plasma is hotter and more highly ionized. Under these conditions, the centroid of the Fe 
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Ka blend and the Fe Kcc/S Hea ratio do not match the observations, even for the lowest values of 
PAM- In Paper II, we suggested that the high degree of mixing in the ejecta that is characteristic 
of all three-dimensional deflagrations is in conflict with the X-ray observations of several Type la 
SNRs, which seem to indicate that Fe and Si are emitting under different physical conditions. This 
conclusion is confirmed by the more detailed work presented here. 

6. MODELING THE SPATIALLY INTEGRATED SPECTRUM 

6.1. Comparing the models to the observed spectrum 

After identifying the most promising explosion models as decribed in § 5, the next step is to 
compare them to the spatially integrated spectrum. In order to perform these comparisons, we 
assume that all the emission from the shocked AM is nonthermal, and describe it with a power 
law continuum of index 2.72 (see discussion in § 2.3). This assumption is based on the results of 
Warren et al. (2005), who used the morphology of the rim emission to constrain the contribution to 
the X-ray flux from a thermal shocked AM component, and found an upper limit of 9% in a subset 
of our region B (226 < 9 < 271). We note that this limit might be somewhat higher in other parts 
of region B where the rim is thicker, but in any case the temperature of the thermal component 
would have to be well below the 2 keV determined by Hwang et al. (1998), because the high energy 
(E>1.5 keV) continuum is known to be dominated by nonthermal emission. Since the exact value 
of the hydrogen column density is unknown, the properties of this low temperature AM component 
would be very difficult to determine from the integrated spectrum. For simplicity, we assume that 
its contribution is negligible. 

As we mentioned in § 3, direct comparison between our models and the observed spectrum is 
not straightforward. Given the limitations of the models and the substantial uncertainties in the 
atomic data included in the spectral codes, it is unrealistic to expect a 'valid fit' from the statistical 
point of view (i.e., with a reduced below the required limit). Under these circumstances, it 
is not appropriate to rely solely on the value of as a measure of the best 'fit', and we use 
additional criteria to decide which models are most satisfactory. In our case, there are only three 
parameters that can be adjusted for each ejecta model: the normalization of the ejecta and AM 
components {norniej, norniAM) and the value of the hydrogen column density Nh- We will adjust 
these parameters in two steps. In a first step, the normalization of the two components will be 
determined by minimizing the statistic in the spectrum between 1.6 and 9.0 keV, with the 
hydrogen column density set to our fiducial value of Nh = 0.6 • 10^^ cm~^. By doing this, we 
guarantee that the high energy spectrum will be approximated to the best ability of the ejecta 
model. In a second step, the normalization parameters will be frozen and the interstellar absorption 
will be adjusted, again by minimizing the statistic, but this time using the spectrum between 
0.8 keV (the location of the brightest line in the Fe L complex) and 10 keV. The goal of this second 
adjustment is to determine whether the ejecta model is capable of reproducing the Fe L/Fe Ka 
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ratio for any given value of Nh- No adjustment whatsoever will be performed on the spectrum 
below 0.8 keV in order to prevent cjccta models with excess emission in this region from forcing 
high values of Nh that might disrupt the Fe L/Fe Ka ratio. 

In addition to the quality of the spectral 'approximation' (as opposed to fit) provided by each 
ejecta model, we employ other consistency checks to evaluate its performance. First, Njj should 
not depart significantly from the values discussed in § 4. Any model requiring Nh to be much 
higher than 0.8 or much lower than 0.4 ■ 10^^ cm~^ is probably overpredicting or underpredicting 
the amount of Fe L emission. Second, the normalization distance Dnorm inferred from the value 
of norniej should be within the 1.5-3.1 kpc range mentioned in § 2.2. The value of Dnorm can be 
determined using the equation 

lOkpc , , 

^norm ~ / ^ V j 

• norniej 

where 10 kpc is the fiducial distance to the source used in the calculation of the ejecta models and 
^ is a correction factor that accounts for the incomplete spatial coverage of region B (which only 
contains photons from ~ 30% of the SNR surface, or 1/^ of the total flux). A more accurate value of 
^ can be found by comparing the fitted line fluxes for the brightest lines in region B listed in Table 
1 to the fluxes for the entire SNR determined by Hwang & Gotthelf (1997). For the four brightest 
line blends, this flux ratio is equal to 2.91 (Si Hea), 2.77 (Si He/3), 2.87 (S He/3), and 2.82 (Ar Hea); 
we adopt a value of 2.8 for ^ as the best approximation to the mean. Third, the value of norniAM 
should make the power law the main contributor to the continuum at high energies, in agreement 
with the results of Warren et al. (2005). Quantitatively, norm am corresponds to the the normalized 
flux at 1 keV from the AM in region B, expressed in units of phot • cm~^ • s~^ • kcV~^. This number 
can be converted to a flux from the shocked AM in the entire remnant using the correction factor ^, 
so that Fam = £,'norniAM should be within the 7.4 — 8. 9 phot • cm~^ • s~^ • keV~^ range determined 
by Fink et al. (1994). Finally, an independent constraint on the distance is provided by matching 
the angular sizes of the fluid discontinuities (RS, CD and FS) determined by Warren et al. (2005). 
Since the dynamics of the FS arc strongly modifled by CR acceleration and the surface of the CD 
is corrugated due to dynamic instabilities, the only straightforward match that we can perform to 
our one-dimensional models with 7 = 5/3 is the angular radius of the RS. We denote by Dug the 
distance required by each model to reproduce the 183" RS radius found by Warren et al. (2005). 



6.2. Results for the spatially integrated spectra 

In Figure 7, we compare the spatially integrated spectrum from region B with eight ejecta 
models selected from our grid. Four of these models have proven the most successful in reproducing 
the diagnostic quantities for the line emission (see § 5.3): DDTa {pam = 2-10~^^ g ■ cm~^, (3 = 0.01), 
DDTc (pam = 2 • 10-24 g . cm-3, /3 = 0.03), 5p0z22.25 (pam = 5 • lO'^^ g • cm-^, (3 = 0.05), and 
PDDa (pam = 2 • lO'^^ g • cm'^, (3 = 0.005). Models DDTc, 5p0z22.25, and PDDa require values 
of P that are off the nodes of the original {pam,P) grid (see Tables 4 and 5). We have selected 
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the values that provide the best possible approximation to the Fe Ka emission (i.e., match the 
observed Fe Ka/S Hea ratio and Fe Ka centroid). Four other ejecta models are also included 
in Figure 7, mainly for illustrative purposes. Model DDTe {pAM = 2 • 10^^"^ g • cm^'^, /? = 0.1) 
completes the sequence of DDT models, and models W7 (pam = 5 • lO'^s g . cm-^, f3 = 0.1), 
SCH {pAM = 5 • 10-25 g . cm-3, (3 = 0.01), and b30_3d_768 {pam = 2 • 10-^5 g . cm-^, /3 = 0.01) 
are provided to assess the viability of one-dimensional deflagrations, sub-Chandrasekhar explosions 
and three-dimensional deflagrations. 

The adjusted values for Nh, norm am, and norruej are listed in Table 6, along with Fam and 
the distance estimates Dnorm and Drs- These two distance estimates are mutually inconsistent 
in all the models presented here. In general, models with higher pam have higher fluxes in the 
shocked ejecta and lower radii of the RS, which result in higher values of Dnorm and lower values 
of Drs, and vice versa (the extreme examples are models 5p0z22.25 and PDDa). We will discuss 
the inconsistency in the distance estimates in § 8.1. 

Delayed detonations (DDTa, DDTc, DDTe, 5p0z22.25) Delayed detonation models are 
clearly superior to the other classes of Type la SN explosions. The low energy X-ray emission 
poses strong constraints on the main parameter in one-dimensional DDT models, the deflagration- 
to-detonation transition density ptr- Model DDTa {pam = 2 • lO-^^g • cm"^, /3 = 0.01), which has 
the highest ptr and therefore synthesizes more Fe and less intermediate mass elements than the 
other models, overpredicts the flux in the Fe L complex, requiring a hydrogen column density of 
0.94 ■ 10^2 cm-2 to compensate this. The shocked Fe is also overionized, as indicated by the excess 
flux coming from the L-shell lines of Fe"'"^'^ and other ions of higher charge around 1.1 keV. Model 
DDTe {pam = 2 • lO^^^ g • cm"^, /3 = 0.1), on the other hand, is probably underpredicting the Fe 
L flux judging from the low Nh that it requires, 0.36 • 10^^ cm-^. Although the shape of the Fe L 
complex is approximately correct in this case, the model clearly overpredicts the O Lya flux at 0.65 
keV, indicating that there is too much O in the outer layers of the ejecta, and that ptr is probably 
too low. The best results are obtained for the intermediate model DDTc {pam = 2 ■ 10"^^ g • cm"^, 
P = 0.03), which can reproduce with reasonable accuracy both the shape of the Fe L complex and 
the flux of the O Lya line for Nh = 0.55 • 10^^ cm-^. The only features that this model does not 
approximate well are the Mg Hea flux at 1.35 kcV and an Fe L-shcll line at 0.7 kcV. A relationship 
between ptr and the properties of the synthetic spectra was already hinted at by Tables 4 and 5, 
because DDT and FDD models with lower ptr required lower values of /3 (see Badenes et al. 2005b, 
for a more general discussion of this topic). The 'off-grid' model 5p0z22.25 {pam = S-IO-^^ g ■ cm"^, 
P = 0.05) is better than DDTa or DDTe, but inferior to DDTc. It requires a hydrogen column 
density of 0.57- 10^^ cm-^, and provides a good approximation to the spatially integrated spectrum, 
fltting the Mg Hea line better than DDTc, but it fails to reproduce the O Lya flux. We note that 
in all the DDT models presented here, the high energy continuum is dominated by the nonthermal 
AM component, with values of Fam between 8.3 and 7.9 • 10"^ phot • cm-^ • s^^ • keV"^, within the 
range determined by Fink et al. (1994). The distance Drs, obtained by matching the location of 
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the RS is within the 1.5-3.1 kpc range in all cases, but the normalization distance Dnorm is always 
too large. 

Pulsating delayed detonations (PDDa) The pulsating delayed detonation model PDDa 
{PAM = 2 • lO^^^g-cm^^, (3 = 0.005) provides a good approximation to the spectrum at high 
energies, but fails at low energies. Although the model requires a reasonable hydrogen column 
density of 0.76 • lO^^cm"^, the prominent Fe+^^ line at ~ 0.8 keV does not appear in the syn- 
thetic spectrum, and neither does the O Lya line. The value of Fam for this model, 7.9 ■ 
10-2 phot • cm 2 ■ s ^ ■ keV ^, is similar to those required by the DDT models, and also within 
the expected range. The low value of pam, however, implies a low emitted flux in the shocked 
ejecta and a large RS radius, which result in an unrealistically low Dnorm 

(0.72 kpc) and a high 

Drs (4.38 kpc). 

One-dimensional deflagrations (W7) In addition to the obvious shortcomings at high ener- 
gies, the W7 {pam = 5 • 10"^^ g • cm"^, /3 = 0.1) model is extremely poor at low energies. This is 
mainly due to the high mass of Mg in the ejecta of W7 (8.5 • IQ-^ Mq, compared to < 10"^ Mq 
in other models). In this case, it is the high flux in the Mg Hea complex, and not an excess of 
Fe L emission, that results in the high value of Nh- The thermal continuum from the W7 model 
is higher than in the DDT and FDD models because of the high amount of unburned C and O in 
the outer layers of ejecta, which leads to a value of Fam below the expected range. The distance 
Dnorm is within the expectations, but Drs is too high. 

Sub-Chandrasekhar Explosions (SCH) The SCH (pam = 5- lO-^^ g • cm-^ /3 = 0.01) model 
is clearly unable to reproduce the observed Fe L/Fe Ka flux ratio. The Fe Ka flux can be increased 
by increasing /3, but this moves the centroid of this line blend to too high energies (see Fig. 6). 
We have verified that an increase in (3 cannot solve problems with Fe L emission either: the flux 
is too high (as indicated by the required Nh of 1.03 • lO^^cm-^) and the overall shape does not 
correspond with the observed spectrum. The AM flux is 7.3 • 10^^ phot • cm-^ • s"^ • keV~^, which 
is just outside the tolerance range derived by Fink et al. (1994). As in the case of model W7, Dnorm 
is within the estimated distance range, but Drs is too high. 

Three-dimensional deflagrations (b30_3d_768) A comparison with the spatially integrated 
spectrum reveals the shortcomings of three dimensional deflagrations even more clearly than the 
diagnostic quantities for the line emission. In the case of model b30_3d_768 (pam = 2-10-^^ g • cm"^, 
f3 = 0.01), the discrepancy between the predicted Fe L/Fe Ka flux ratio and the observed spectrum 
is dramatic, with the model requiring a hydrogen column density of 2.33 • 10^^ cm-^. In fact, it 
can be said that the spectral adjustment is poor everywhere, and it is clear from Figure 6 that no 
combination of pam and (3 is going to improve it. The high level of thermal continuum in the ejecta 
model, which comes mostly from the large amounts of unburned C and O in the outer ejecta layers. 



-20- 



forces a very low value of Fam (1-0 ■ 10~^ phot • cm~^ ■ ■ keV"-*^) that is impossible to reconcile 
with the Ginga observations. As in most models that require a low value of pam, Dnorm is too low 
and Drs is too high. 

To conclude this section, we mention an unrelated, but interesting fact about the high energy 
continuum that we have noticed during our attempts to approximate the integrated spectrum of 
the Tycho SNR. Hwang et al. (2002) and Warren et al. (2005) found that the featureless rim emis- 
sion extracted from Chandra observations could be fitted equally well by thermal and nonthermal 
models. In our integrated XMM-Newton spectra, however, we have found that thermal models for 
the continuum always underpredict the flux above 8 keV (e.g., see Figure 2 in Bravo et al. 2005). 
This is revealed by the higher effective area of XMM at high energies, and it constitutes yet another 
argument in favor of the predominantly nonthermal origin for the shocked AM emission advocated 
in Warren et al. (2005). 

7. SPATIAL MORPHOLOGY OF THE LINE EMISSION 

Although the focus of this paper is on the spatially integrated X-ray emission from the Tycho 
SNR, the spatial distribution of this X-ray emission also contains interesting information that can 
shed light on the properties of the SN ejecta and the SNR shocks. In Paper II, for instance, 
we proposed that partial collisionlcss heating at the reverse shock could explain why the Fe Ka 
emission peaks interior to Fe L and Si Hea in both Tycho (Hwang & Gotthelf 1997) and Kepler 
(Cassam-Chenai et al. 2004). We note that partial collisionlcss electron heating has been required 
by all the models that have shown some level of success in reproducing the spatially integrated 
spectrum of the Tycho SNR in the previous sections. We want to conclude our study by doing a 
preliminary comparison between our models and the spatial distribution of the line emission in the 
Tycho SNR. 

We have chosen two models for this comparison: DDTc {pAM = 2 ■ 10~^^g • cm~^, /? = 0.03) 
and b30_3d_768 {pam = 2 • 10~^^ g • cm~^, /? = 0.01). Both models have a similar amount of partial 
collisionless electron heating at the RS, but their ejecta structures are very different (layered vs. 
well-mixed) , and the order of magnitude difference in pAM places them at very different evolutionary 
stages. At the age of Tycho, the RS has heated ~ 1.0 Mq of ejecta in model DDTc, but only 
^ 0.5 M0 in model b30_3d_768. However, physical conditions in the shocked ejecta are qualitatively 
similar (see Figs. 8 and 9): p and Ugt increase monotonically from RS to CD, but Tg has local 
maxima at the CD (due to coUisional heating in the plasma over time) and close to the RS (due to 
collisionless heating at the RS). 

Taking advantage of the superior spatial resolution of Chandra, we have produced radial profiles 

for the emission in three spectral bands, using data from the ACIS-I CCD detectors: 0.8-0.95 keV 
(Fe L), 1.63-2.26 keV (Si K, including Si Hea, Si He/3, and Si Lya), and 6.10-6.80 keV (Fe Ka). 
The data are from the ~ 150 ks observation described in Warren et al. (2005), and were extracted 
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from the same region B whose spatially integrated spectrum we have used throughout the present 
work. For the Fe Ka profile, the underlying continuum was subtracted as in Warren et al. (2005). 
In Figure 10, we plot the three line emission profiles alongside the predictions from the two models. 
The model profiles have been scaled so the position of the RS coincides with the 183" radius found 
by Warren et al. (2005) (first vertical dotted line in Fig. 10), which is equivalent to placing each 
model at the distance Drs listed in Table 6. The model profiles have also been normalized to the 
peak value in the data for each energy band. Prom the spectra presented in Fig. 7, it is clear that 
model DDTc produces the correct spatially integrated flux in each of the three spectral bands (for 
Nh = 0.55 • 10^2 cm-2), but model b30_3d_768 does not. 

Our one-dimensional models look very different from the data, in part because they do not 
extend to the observed radius of the CD, with the radius of the outermost ejecta falling ~ 10% 
short for model DDTc and ~ 20% short for model b30_3d_768. There are several processes that 
we have not included in our simulations and could increase the radial extent of the ejecta, but the 
morphology of the CD and its power spectrum seem to indicate that Rayleigh- Taylor instabilities 
(RTI) are the dominant mechanism in the Tycho SNR (see § 4, § 7.3 and Fig. 6 in Warren et al. 
2005). Wang &: Chevalier (2001) explored the development of RTI in the context of an exponential 
density profile (which is considered the best generic approximation to Type la SN ejecta, see 
Dwarkadas & Chevalier 1998) interacting with a uniform AM, and found an increase of ~ 10% in 
the radius of the outermost ejecta at the age of the Tycho SNR (see their Figure 4). This would be 
sufficient for model DDTc, but not for model b30_3d_768. Despite the obvious limitations of our 
one-dimensional models, wc note that model DDTc gives roughly the right morphology: the Fc Ka 
profile peaks interior to the other two, the maxima of the Fe L and Si K profiles are at the right 
locations, and the Si K-shell emission extends a little bit beyond the Fe L-shcU emission. This is 
the result of the combination of the layered structure of the ejecta with the distribution of physical 
conditions in the shocked plasma represented in Fig 8. The peaks in the Fe L and Si K emission 
coincide with the most dense layers with high Fe and Si abundances. The Fe Ka emission, on the 
other hand, has a strong temperature dependence, and is greatly enhanced in the hot region close 
to the RS produced by the collisionless electron heating at the shock front. For model b30_3d_768, 
on the other hand, all the predicted profiles have the same morphology. In this case, the peak of 
the Fe Ka emission coincides with the other two because the outermost ejecta layers are rich in Fe. 
A higher amount of collisionless electron heating will enhance the Fe Ka emission close to the RS, 
but this would not be beneficial for the spatially integrated spectrum of this model (see Fig. 6). 



8. DISCUSSION 

In the present work we have compared the X-ray emission from the shocked ejecta in the Tycho 
SNR to the predictions from our models based on one-dimensional hydrodynamic simulations and 
NEI calculations. We have focused on the spatially integrated X-ray spectrum, but wc have also 
taken into consideration the spatial morphology of the line emission and the dynamics of the 
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underlying hydrodynamic calculations. We have shown that our best models achieve a remarkable 

level of success in explaining the global properties of the ejecta emission. The only inconsistency 
that wc have found in our comparisons is the conflict between the distance estimates Dnorm and 
Djis- In this section, we will discuss this conflict and the impact that it has on our results, and 
we will summarize the case for a delayed detonation as the physical mechanism responsible for SN 
1572. 

8.1. The distance estimates: Dnorm vs. Drs 

The distance estimates Dnorm and Drs never agree with each other in the models listed in 
Table 6 or, for that matter, in any of the models that wc have tested. In the most successful class 
of models (the delayed detonations), Djis is always within the expected 1.5-3.1 kpc range, while 
Dnorm IS Systematically larger, with values ranging between 3.58 and 6.95 kpc. Some insight into 
the origin of this discrepancy can be gained by inspecting the radial line brightness profiles given in 
Figure 10. The volume of X-ray emitting ejecta can be estimated from the measured locations of the 
RS and CD given by Warren et al. (2005) (vertical dotted lines in Fig. 10). This volume is always 
considerably larger than the volume occupied by the ejecta in our one-dimensional models, assuming 
the RS radius in the models matches that of Tycho. As wc have discussed in the previous section, 
this is probably due to the effect of Ray Icigh- Taylor instabilities at the CD. Whatever the cause, 
spreading a fixed mass of shocked ejecta over a larger volume will decrease the average density, 
reducing the intrinsic emissivity of the model and lowering the inferred normalization distance 
Dnorm- The value of DrSi being determined by the hydrodynamic evolution, should remain nearly 
the same (see § 3.1 in Wang & Chevalier 2001). 

The specific details of how the development of RTI will affect the emitted flux, the properties of 
the spatially integrated spectrum and the spatial morphology of the line emission are complex, and a 
detailed study of these issues is outside the scope of this paper. However, we want to point out here 
that the difi"ercnces in the spatially integrated synthetic spectrum between a one-dimensional model 
and a multi-dimensional model including RTI might be relatively small. The spatially integrated 
emission is determined by the spectrum emitted by each individual fluid element in the SN ejecta, 
which only depends on its composition, its preshock density (or equivalently, the time at which it is 
overrun by the RS) and its postshock evolution. The dynamics of the RS do not change significantly 
in multi-D simulations including RTI (see Chevalier et al. 1992; Wang &: Chevalier 2001), so the 
density at which a fiuid clement with a given composition is shocked will not be altered. After 
the passage of the RS, the development of RT fingers will enhance the adiabatic expansion of the 
fluid element, reducing the density in the postshock evolution. The extent of this reduction can be 
estimated by looking at Figure 4 in Wang &; Chevalier (2001): at the age of the Tycho SNR, the 
angle-averaged density profile affected by RTI only shows a factor 2 decrease with respect to the 
highest density ejecta in a one-dimensional simulation (which is just behind the CD). In this region, 
the high density RT fingers penetrate into the low density shocked AM, with the latter component 
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filling most of the volume (see Figure 2 in Wang &: Chevalier 2001), so the actual drop in density 
of the shocked ejecta should be less than a factor 2. We know that the emissivity and emission 
measure of an X-ray plasma scale with the square of the density, while the coUisional processes in 
the plasma (ionization, electron heating) only have a linear density dependence (see Eqns. 1, 2, 
and 6 in paper I). Therefore, it does not seem unreasonable to think that the flux emitted by the 
shocked ejecta may drop by ~ 50% (driving Dnorm in the DDT models into the expected range), 
without major changes being introduced to the value of ngt and Tg inside each fluid clement and 
hence to the properties of the spatially integrated spectrum. In this context, the differences in the 
synthetic spectra from one Type la SN explosion model to another should be conserved to a large 
extent. However, we insist that detailed multi-dimensional simulations are required to investigate 
these issues. We plan to perform such simulations in a forthcoming paper. 

We conclude that the discrepancy between Dnorm and Drs that we have encountered might be 
an intrinsic feature of one-dimensional calculations. Prom what we know about the structure and 
growth of RT fingers at the CD, their effect on the properties of the synthetic spectra in terms of 
ionization timescalcs and electron temperatures might be limited. Prom a practical point of view, 
the presence of RTI at the CD might modify the optimum value of pAM for a given Type la SN 
explosion model, but it is very unlikely that it will make viable the synthetic spectrum of a model 
that did not prove valid in a one-dimensional calculation. 

8.2. A delayed detonation as the origin of SN 1572 

The main result of our work is that one-dimensional delayed detonation models are successful 
in reproducing the fundamental properties of the X-ray emission from the shocked ejecta in the 
Tycho SNR. None of the other explosion paradigms that we have tested (prompt detonations, 
sub-Chandrasckhar explosions, pulsating delayed detonations, one-dimensional deflagrations and 
three-dimensional deflagrations) shows a comparable level of success, and most of them fail even 
simple comparisons to the observations. The delayed detonation models reproduce 12 out of the 
13 diagnostic quantities for the line emission (see Figs. 4, 5, and 6; and Tables 4 and 5), and they 
provide good approximations to the spatially integrated XMM-Newton spectrum at high energies 
(see fig.7). All the delayed detonation models that we have tested require similar values of pam 
(2 • 10~^^ g • cm~^ for our grid DDT models, 5 • 10~^^ g • cm~^ for the 'off-grid' model 5p0z22.25) 
and P (between 0.01 and 0.1). These values lead to estimates for the distance Djis, the nonthermal 
flux from the shocked AM Fam, and (in the case of model DDTc) the hydrogen column density 
Nh, which are fully consistent with the available observational constraints. We emphasize that, 
by taking these estimates into account, we are evaluating our models from a global point of view, 
not based on the properties of their synthetic spectra alone. We have also shown that the spatial 
morphology of the line emission predicted by delayed detonation models with partial collisionless 
electron heating is roughly consistent with the Chandra observations of the Tycho SNR, at least 
within the limitations of one-dimensional models. The inconsistency between the distance estimates 
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from matching the RS radius (Drs) and from the normahzation of the spectrum (Dnorm) that we 
found in all the models could also be related to the limitations of our one-dimensional hydrodynamic 
calculations. A detailed study of the effect of the Rayleigh- Taylor instability on the X-ray emission 
from the shocked ejecta would be necessary to gain further insight on the spatial morphology of 
the line emission and the values of Dnorm- 

We believe that the unavoidable limitations of one-dimensional simulations do not compromise 
our results concerning delayed detonation models (see discussion in 8.1). This is specially true 
because we have found no alternative paradigm with an acceptable level of success whose predictions 
could be reconciled with the observations by invoking the limited modifications that we expect 
from multi-dimensional effects like RTI. We note that our conclusions do not agree with those of 
Sorokina et al. (2004), who show a preference for three-dimensional deflagrations for the Tycho 
SNR. However, these authors only consider two deflagration models (W7 and a model similar to 
b30_3d_768), and their results are based on a very simple qualitative comparison between models 
and observations. To evaluate the merits of our slightly different approaches on an equal footing, 
synthetic spectra from delayed detonation models should be produced using the methods of Sorokina 
et al., and they should be compared to the observations in a rigorous, quantitative way. 



8.3. The case for model DDTc {pam = 2 ■ 10~^^ g • cm^^, (3 = 0.03) 

Among the delayed detonation models that we have tested, model DDTc {pAM = 2-10~^^ g • cm" 
[3 = 0.03) gives by far the best approximation to the spatially integrated X-ray spectrum. The 
performance of this model is surprisingly good: despite the limitations of our one-dimensional cal- 
culations and the large uncertainties in the atomic data used to generate the synthetic spectra, the 
model can reproduce the fluxes and centroids of almost all the lines from O, Si, S, Ar, Ca and Fe 
in the shocked ejecta between 0.6 and 7.0 kcV (see Fig. 7). Our preference for this model is based 
on the constraints posed by the spectrum at low energies (below 1 keV) on the main parameter 
involved in delayed detonation explosions, the transition density from deflagration to detonation 
Ptr- Model DDTc matches both the Fe L/Fe Kct flux ratio and the O Lya emission for a reasonable 
value of Nh (0.55 • 10^^ cm~^), unlike models with higher or lower pt^ such as DDTa or DDTe. The 
second most successful delayed detonation model, 5p0z22.25, although calculated using a different 
SN explosion code, is very similar to DDTc. Both models have similar values of ptr (2.2-10^ g • cm~^ 
for DDTc and 2.5 • 10^ g • cm^'^ for 5p0z22.25), similar kinetic energies (1.16 and 1.20 • 10^^ erg), 
and they synthesize similar amounts of most elements (Fe: 0.80 and 0.74 Mq, Si: 0.17 and 0.22 
M©, S: 0.13 and 0.12 Mq). 

The fact that the model agrees so well with the observed spectrum implies that the bulk 
properties of the synthetic spectrum (i.e., emission measures and emission measure averaged tem- 
peratures and ionization timescales) are roughly correct for the most prominent elements, and can 
be considered representative of the state of the shocked ejecta in the Tycho SNR. The spectral 
adjustment provided by the model, however, is not perfect. As we discussed in § 5.2, the centroid 
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of the Ca Hea blend is underpredicted by the model. A close examination of the residuals in Fig. 
7 reveals that the centroids of the Si Hea, S Hea, Ar Hea and Fe Ka blends in the synthetic spec- 
trum are also offset towards lower energies with respect to the data. The shape of these residuals 
is suggestive of a gain shift, but we have verified that this is not the case - the shift is indeed real, 
and indicates that the plasma in the model is slightly underionized. Increasing the value of pam 
gradually between 2.0 and 5.0 • 10~^^ g • cm~^ might solve this issue, but given the limitations of 
one-dimensional models discussed in § 8.1, we do not believe that such refinements are justified. 

The mean velocity of the shocked ejecta in this model is ~ 2 • 10^cm-s~^. For a radius 
of ~ 7 • 10^^ cm (see Fig. 8), this gives an expansion of ~ 0.09% yr~^. This is in very good 
agreement with the expansion parameters of the interior features of the SNR measured in radio 
{rj ~ 0.44, Reynoso et al. 1999) and X-rays {rj ~ 0.45, Hughes 2000), which translate into a 
percentual expansion of of ~ 0.11% yr~^. We note that, when the peaks of the Si Hea blend are 
aligned, the width of the blend seems larger in the data than in the synthetic spectrum, suggesting a 
Doppler broadening of the line emission in the shocked ejecta. This would certainly have interesting 
implications for the dynamics of the Tycho SNR, but a quantitative analysis of this effect would be 
challenging given the spectral resolution of the CCD cameras onboard Chandra and XMM-Newton. 

We conclude our comments with a reference to the historical observations of SN 1572. The 

light curve calculated for model DDTc yields a bolometric magnitude at light curve maximum 
M max of -19.51 mag (see Table 3 in Paper II), which is at the upper limit of the M-v — —19.24 — 
5/o(/(I?/3.0kpc) ± 0.42 range derived by Ruiz-Lapuente (2004) taking the value of 2.8 kpc for D 
proposed by this author, but 0.17 mag too bright at Drs = 2.59 kpc. Given the fact that the 
estimates of Ruiz-Lapuente are based on naked eye observations, and the large range of peak 
bolometric magnitudes in our Type la SN model grid (more than 2 mag, see Table 1 in Paper I), 
this can be considered a good agreement. An explosion that synthesizes ~ O,8M0 of ^^Ni (like 
model DDTc) is consistent with the hypothesis that SN 1572 was a 'normal' Type la SN, and not 
a subluminous event. 



9. CONCLUSIONS 

In the present work, we have compared synthetic X-ray spectra for the ejecta emission in Type 
la SNRs with the X-ray observations of the Tycho SNR, with the goal of constraining the explosion 
mechanism that gave birth to the supernova of 1572. In our comparisons, we have found a very 
strong bias towards delayed detonation models, to the point that no other Type la SN explosion 
paradigm has even come close to reproducing the fundamental properties of the X-ray emission 
from the ejecta in the Tycho SNR. This result has a high significance, and we believe that it is not 
affected by the limitations of our one-dimensional models. 

Among the delayed detonation models that we have tested, we obtained the best results for 
model DDTc, expanding into a uniform AM of pam = 2 • 10~^^ g • cm^"^, and with a modest amount 
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of collisionless heating at the RS (/? = 0.03). At the age of the Tycho SNR, the synthetic spectrum 
calculated for this model can account for almost all the fundamental properties of the ejecta emission 
in the XMM-Newton observation, including the line emission from O, Si, S, Ar, Ca, and Fe between 
0.6 and 7.0 keV. In agreement with the results of Warren et al. (2005), the spectrum only requires 
a nonthermal contribution from the AM, with a flux of 8.1 ■ 10~^ phot • cm"^ • • keV~^, which is 
within the range determined by Fink et al. (1994), and a hydrogen column density of 0.55-10^^ cm~^, 
which is also compatible with previous observations (e.g. Hwang et al. 1998). By matching the 
radius of the RS in the model to the 183" radius found by Warren et al. (2005), we obtain a 
distance estimate of 2.59 kpc, within the range determined by Smith et al. (1991). The distance 
estimated from the normalization of the emitted spectrum is not consistent with this, but the origin 
of this discrepancy is probably in the inherent limitations of the one-dimensional hydrodynamic 
calculations that underlie our synthetic spectra. The spatial morphology of the line emission 
is qualitatively consistent with the Chandra observations of the Tycho SNR, although multi-D 
simulations are required for a quantitative comparison with observations. We conclude that model 
DDTc represents the closest approximation to the structure of the ejecta in the Tycho SNR that 
we can find within our grid of synthetic spectra. The deflagration-to-detonation transition density 
of this model is ptr = 2.2 • 10^ g ■ cm~^, its kinetic energy is = 1.16 • 10^^ erg, and it synthesizes 
0.12 M© of C and O, 0.17 M© of Si, 0.13 M© of S, 0.033 M© of Ar, 0.038 M© of Ca, and 0.80 M© 
of Fe. We note that a future analysis based on more sophisticated (e.g., multi-D) SNR models 
and/or a different grid of Type la explosion models might lead to a revision of the specific values 
for the kinetic energy and nucleosynthetic yields that we cite here, but we do not anticipate major 
changes. 

Our results have interesting implications for the analysis of X-ray observations of Type la 
SNRs. Our models allow to consider the emission from all the elements in the ejecta at once, 
and they are characterized by very few parameters with a clear physical meaning: the age of 
the SNR t, the AM density pam, the efficiency of collisionless electron heating at the RS (3, and 
the explosion mechanism. These models constitute a promising way to approach the analysis of 
the ejecta emission, foregoing many of the ambiguities inherent to more conventional techniques 
based on spectral fitting (for a discussion, see Rakowski et al., ApJ, submitted). Futhermore, by 
interpreting the X-ray spectrum of the SNR in terms of realistic supernova explosion models, it is 
possible to distinguish between remnants originated by core collapse or thermonuclear supernovae 
with a higher level of confidence. In future papers, we plan to apply our grid of synthetic spectra 
to the other prototype Galactic Type la SNR, SN 1006, and to improve our models by including 
multi-dimensional effects like Rayleigh-Taylor instabilities. 

The implications for the study of SN explosions arc also interesting. The X-ray observations 
of SNRs open a new window onto the structure of the SN ejecta, totally independent of the optical, 
ultraviolet and infrared light curves and spectra of SNc. The two approaches are complementary, 
because they involve completely different challenges and constraints: modeling the radiative transfer 
in the expanding ejecta in the case of the spectral evolution of Type la SNe (see Branch et al. 2005, 
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and references therein) vs. modeling the dynamics of the shocked plasma in Type la SNRs. One 
advantage of SNRs is that the emitting plasma is optically thin, and all the elements can be seen 
in the X-rays at the same time. This makes the success of delayed detonation models for the Tycho 
SNR all the more significant, and the agreement with previous studies of Type la SNe, which also 
showed a preference for delayed detonations (e.g. Hoflich & Khokhlov 1996), more compelling. It is 
worth noting that we have only examined the properties of the ejecta emission from the Tycho SNR 
on the western side of the remnant, where they are reasonably homogeneous. Significant spatial 
inhomogeneities exist in the eastern side (Vancura et al. 1995), which could provide important clues 
about asymmetries in Type la SN explosions. 

We conclude with a reminder that, despite their apparent success, delayed detonations are 
just a phenomenological model for Type la SNe. Since this explosion paradigm was first proposed 
by Khokhlov (1991) ^, all attempts to provide a physically consistent mechanism that can explain 
the transition from deflagration to detonation have failed. Some simulations have explored delayed 
detonation explosions in three dimensions (Gamezo et al. 2005) , but without providing any further 
insight into the physical mechanism that governs the transition. Yet, our results provide unambigu- 
ous evidence that the ejecta of the Type la SN that originated the Tycho SNR are stratified, which 
implies that a detonation was involved at some stage in the physics of the explosion. Very recently, 
two new explosion mechanisms for Type la SNe have been proposed that lead to detonations in a 
self-consistent way: gravitationally confined detonations (Plewa et al. 2004) and pulsating reverse 
detonations (Bravo & Garcia-Senz, ApJ, submitted). Both these paradigms result in a stratified 
structure for the SN ejecta, which is qualitatively compatible with our results, but a quantitative 
comparison is required to evaluate their merits in the context of the X-ray emission from Type la 
SNRs. There is clearly much work left to be done before the physics of Type la SN explosions can 
be completely understood, and we expect the X-ray observations of young, ejecta-dominated Type 
la SNRs to make a significant contribution to this crucial issue. 

This work is dedicated to the memory of Vicente Badenes, the grandfather of Carles, who 
passed away on March 9, 2005. 
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Fig. 1.— The XMM-Newton EPIC MOS image of the Tycho SNR, with the spatial extraction 
regions A and B indicated. 
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Fig. 2. — XMM-Newton EPIC MOSl spectra of regions A and B. The spectrum of region B has 
been offset for clarity, dividing the flux by 5.0. The most important lines and line complexes have 
been marked. 
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Fig. 3 — Fits to the line emission between 1.65 and 9 keV in the XMM-Newton EPIC MOSl 
spectrum of regions A and B. The most important lines and line blends have been labeled for 
clarity. 
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Fig. 4. — Comparison between the predicted line flux ratios and centroids and the observed values 
for models DDTc (left) and PDDa (right). In each case, the thirteen panels correspond to the eight 
diagnostic line flux ratios (top two rows) and the five diagnostic line centroids (bottom row). The 
predicted values arc represented in the horizontal axis as a function of pAM, which is mapped on 
the vertical axis. The values of /? are plotted with different styles: solid line and diamonds for 
Pmin, dotted line and triangles for /? = 0.01, and dashed line and squares for /? = 0.1. The observed 
values are represented by vertical dash-triple-dotted lines, with the borders of the tolerance ranges 
described in § 5.2 as dotted lines and the 'allowed' regions shaded in gray. Absent points denote 
lines that are too weak to derive a flux or a centroid. 
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Fig. 5. — Comparison between the predicted line fiux ratios and centroids and the observed values 
for models DEFc (left) and W7 (right). See Figure 4 for an explanation of the plots and labels. 
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Fig. 6. — Comparison bet"ween the predicted line flux ratios and centroids and the observed values 
for models SCH (left) and b30_3d_768 (right). See Figure 4 for an explanation of the plots and 
labels. 
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Fig. 7. — Comparison between the ejecta models and the spatially integrated spectrum of region 
B. The ejecta model, with the corresponding values pam (in units of 10~^^ g ■ cm~^) and /3, is 
indicated in each plot, as well as the adjusted value of Nh- The power law component is displayed 
for clarity (dotted line). 
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Fig. 8. — The radial structure of tire sliocked ejeeta in model DDTc {pAM = 2 • 10~^^ g • cm~'^, (3 = 
0.03) at the age of the Tycho SNR. (a) Radial distribution of density p and specific internal energy 
e; (b) mean number of electrons per ion, Z, with an indication of the ejecta layers dominated by 
Fe, Si-S, and C-0; (c) electron and ion temperatures; (d) ionization timescale. The positions of the 
RS and CD are outlined by the limits of the temperature plots in (c) . 




Fig. 9.— Same as Fig.8, but for model b30_3d_768 {pam = 2 • 10"^^ g • cm^^ p = 0.01). 
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Fig. 10. — Radial line emissivity profiles for Fe L, Si K, and Fe Ka in models DDTc (left) and 
b30_3d_768 (right), compared to the Chandra profiles for region B. The value of pam in the models 
is given in units of 10~^^ g • cm~^. The locations of the RS and CD (from Warren et al. 2005) have 
been indicated by dotted lines. 



Table 1. Line fluxes and centroids in the Tycho SNR 











Region A 




Region B 






Expected 


Centroid 


Flux 


Centroid 




Flux 


Line 


Ion and Transition 


Centroid (keV) 


(keV) 


(10~^phot ■ cm~2 ■ s~l) 


(keV) 


(10- 


"^phot • cm~^ • s~^) 


Si Hea 


Si+12^ n = 2^n=l 


~ 1.86 


l.»b»4_(, QQQ2 


30. 6+01 


1.00(0_Q QQQ2 




18.0 ±0.1 


Si He/3 


Si+12, ls3p ^ ls2 


2.182 






9 1S4+0002 

^.i»4_(, 




1.54 ± 0.04 


Si Hc7 


Si+12, ls4p ^ 


2.294 


b 


0.35 X SiHc/3 






0.35 X Si He/3 


Si Lya 


Si+13, 2p ^ Is 


2.006 




< 1.04 






< 0.58 


Si Ly/3 


Si+13, Is 


2.377 




0.14 X Si Lya 






0.14 X SiLyo 


S Hea 


S+14, n = 2^n = l 


~ 2.45 


2.447 ± 0.001 




2.447 ± 0.001 




0-0.08 


S He/3 


S+", ls3p^ ls2 


2.884 




0.52 ± 0.03 






0.29 ±0.02 


S Hc7 


S+", ls4p ^ ls2 


3.033 




0.40 X S Hc/3 






0.40 X SHe/3 


S Lya 


2p ^ Is 


2.623 




< 0.05 






< 0.02 " 


Ar Hea 


Ar+"\ n = 2 ^ n = 1 


~ 3.1 


n ,01+0.003 

0.101_Q 


n fis+O-04 

U-O8_o.03 


o 19^^+0.003 




0.38_o ,j2 


Ar He/3 


Ar+ifS, ls3p ^ Is^ 


3.685 




0.06 it 0.02 






02+° "^ 


Ar He7 


Ar+16, is4p -* ls2 


3.875 




0.57 X ArHe/3 






0.57 X ArHe/3 


Ca Hea 


Ca+i*, n = 2 ^ n = 1 


~ 3.88 


3.879 ± 0.009 


0.24 ± 0.03 


o SK7+0 013 




0.15 ±0.02 


Fe Ka 


Several, n = 2 — > n = 1 


~ 6.45 


6.462 ± 0.009 


0.26 ±0.02 


6.456 ±0.013 




0.15 ±0.02 



^The limits given are the formal 90% confidence ranges (Ax^ = 2.706). 

''Centroids marked as ... were not fitted. 

'^For tiie Si and S Lya lines, only the 3<7 upper limits to the flux are given. 
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Table 2. Diagnostic line ratios for the Tycho SNR 



Line Ratio 


Fitted Values 




Region A 


Region B 


SiHe(/3 + 7)/SiHea 


0.115 ±0.003 


0.116 ±0.004 


SiLya/SiHea 


< 0.034 


< 0.032 


SHe(/3 + 7)/SHea 


0.090 ± 0.006 


0.084 ± 0.008 


SLya/SHea 


< 0.006 


< 0.004 


SiHea/SHeo 


3.76 ± 0.07 


3.80 ±0.09 


ArHeo/SHea 


0.084 ± 0.006 


0.080 ± 0.007 


CaHca/SHca 


0.029 ± 0.004 


0.032 ±0.006 


FeKa/SHea 


0.032 ± 0.003 


0.032 ± 0.004 



■^For line fluxes with asymmetric 90% confidence 
ranges, the error in the line flux ratios has been cal- 
culated assuming symmetric confidence ranges with the 
largest of the two deviations. 



Table 3. Extraction energy windows for the model spectra 



Line or blend 


Energy Window (keV) 


Si He« 


1.80 - 1.90 


Si He/3 


2.17-2.20 


Si He7 


2.28 - 2.32 


Si Lya 


1.97 - 2.02 


S Hea 


2.35 - 2.50 


S He/3 


2.87-2.90 


S He7 


3.03 - 3.05 


S Lya 


2.61 - 2.64 


Ar Hea 


3.05-3.16 


Ar He/3* 


3.67 - 3.70 


Ca Hea'' 


3.70 - 3.93 


Fe Ka" 


6.20 - 6.90 



''Might include some Ca Ka lines around 
3.7 keV. 

''Might include some Ca Ka lines around 
3.7 keV and a contribution from Ar He7 at 
3.88 keV. 



Might include some Fe Hea lines 
around 6.67 keV. 



Table 4. Models vs. observations: line flux ratios 



Parameters Line Flux Ratios 





PAM 




Si He(/3 + 7)/ 


Si Lya/ 


SHe(/J + 7)/ 


SLya/ 


Si Hea/ 


Ar Hea/ 


Ca Hea/ 


Fe Ka/ 


Model 


(10-24 g . cm-3) 




Si Hea 


Si Hea 


S Hea 


S Hea 


S Hea 


S Hea 


S Hea 


S Hea 


DDTa 


2.0 


0.01 






/ 


/ 






/ 


/ 


DDTc 


2.0 


0.01 -f- 0.1 


/ 




/ 


/ 






/ 


/ 


DDTe 


2.0 


0.1 






/ 


/ 






/ 


/ 


5p0z22.25 


5.0 


0.01 ^ 0.1 




■/ 


/ 


/ 




■/ 






PDDa 


0.2 


/3min ^ 0.01 




^ 




/ 




/ 


/ 


/ 


PDDc 


0.2 


0.01 






/ 


/ 




-/ 


■/ 




PDDe 


0.2 


0.1 




■/ 




t 


■/ 








DEFa 


5.0 


0.1 




t 




t 








/ 


DEFc 


5.0 


0.01 




t 




t 








/ 


DEFf 


2.0 


0.01 ^ 0.1 




t 




t 










W7 


0.5 


0.1 






/ 


/ 








/ 


SCH 


0.5 


0.01 






/ 


/ 




t 






DET 


2.0 






i 






$ a 






$ a 


B30U 


2.0 


0.1 




■/ 






t 


i 




t 


b30.3d.768 


0.2 


0.01 




■/ 








t 




■/ 


Note. — 


within tolerance 


range; ~; marg 


;inal (on the border of the tolerance range); f: 


discrepant; ^: line is 


too weak to calculate 


a flux or is 



absent from the synthetic spectrum. 

''The Si Hea and Fe Ka blends are both present in the synthetic spectrum, but the S Hea blend is too weak to calculate these line flux ratios. 



Table 5. Models vs. observations: line centroids 





Parameters 






Line Centroids 




Obscrvables 




PAM 


P 












Jricproduced 


Model 


(10-24 g.cm-3) 




Si Hea 


S Hea 


Ar Hea 


Ca Hea 


Fe Ka 


(out of 13) 


DDTa 


2.0 


0.01 




/ 


■/ 






12 


DDTc 


2.0 


0.01 0.1 






■/ 


t 




12 


DDTe 


2.0 


0.1 


/ 




■r 


t 


t 


11 


5p0z22.25 


5.0 


0.01 ^ 0.1 






•r 






12 


PDDa 


0.2 


Pmin ^ 0.01 




/ 


•/ 




/ 


12 


PDDc 


0.2 


0.01 


/ 




/ 


t 


t 


11 


I'DDe 


0.2 


0.1 


/ 




/ 


t 


t 


y 


DEFa 


5.0 


0.1 






•/ 


t 




9 


DEFc 


5.0 


0.01 






■/ 




t 


9 


DEFf 


2.0 


0.01 -r 0.1 




/ 




t 




9 


W7 


0.5 


0.1 




/ 




t 


t 


9 


SCH 


0.5 


0.01 




/ 




t 




8 


DET 


2.0 


Pmin 




$ 


t 


t 




2 


B30U 


2.0 


0.1 






i 


i 


t 


6 


b30.3d.768 


0.2 


0.01 






i 


i 


t 


8 



Note. — ■/ : within tolerance range; ~; marginal (on the border of the tolerance range); f: discrepant; ^: line is 
too weak to calculate a centroid or is absent from the synthetic spectrum. 
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Table 6. Type la explosion models for the Tycho SNR 











normAM 




Fam 






Model 


PAM 




Nh " 


(IQ-^phot • cm^-^- 


(phot ■ cm— 


(10-2phot • cm"2. 


Dnorm 


Drs 




(10-24 g • cm-3) 




(10^2 cm-2) 


-keV^i) 


s-l -keV-l) 


s-i -keV-i) 


(kpc) 


(kpc) 


DDTa 


2.0 


0.01 


0.94 


2.97 


2.79 


8.3 


3.58 


2.57 


DDTc 


2.0 


0.03 


0.55 


2.89 


1.78 


8.1 


4.48 


2.59 


DDTe 


2.0 


0.1 


0.36 


2.81 


1.36 


7.9 


5.12 


2.54 


5p0z22.25 


5.0 


0.05 


0.57 


2.89 


0.74 


8.1 


6.95 


1.95 


PDDa 


0.2 


0.005 


0.76 


2.81 


68.3 


7.9 


0.72 


4.38 


W7 


0.5 


0.1 


1.07 


2.15 


6.37 


6.0 


2.36 


3.54 


SCH 


0.5 


0.01 


1.03 


2.60 


10.82 


7.3 


1.82 


3.44 


b30.3d.768 


0.2 


0.01 


2.33 


0.35 


59.8 


1.0 


0.77 


3.41 



^The spectral approximations presented in Fig. 7 and in this Table are not statistically valid fits (see text), which makes it hard to 
define confidence ranges for the adjusted parameters Njj, normAM, and norniej- We do not provide any estimates on the errors for 
these parameters, but note that they may be large. 



